Descritpion of Exotic Nuclei Using Continuum Shell Model 
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(today) 

In weakly bound exotic nuclei, number of excited bound states or narrow resonances is small and, 
moreover, they couple strongly to the particle continuum. Hence, these systems should be described 
in the quantum open system formalism which does not artificially separate the subspaces of (quasi- 
) bound and scattering states. The Shell Model Embedded in the Continuum provides a novel 
approach which solves this problem. Examples of application in sd shell nuclei will be presented. 
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I. INTRODUCTION 

A realistic account of the low-lying states properties in 
exotic nuclei requires taking into account the coupling be- 
tween discrete and continuum states. This aspect is par- 
ticularly important in studies near the drip line where one 
has to use both structure and reaction data to understand 
basic properties of these nuclei. Within the Shell Model 
Embedded in the Continuum (SMEC) [0, one may ob- 
tain a unified description of the divergent characteristics, 
such as the spectra (energies of states, transition proba- 
bilities, proton/neutron emission widths, (3— decays, etc.) 
and the reactions involving one-nucleon in the contin- 
uum (proton/neutron capture processes, Coulomb disso- 
ciation reactions, elastic/inelastic proton/neutron reac- 
tions, etc.). The accumulation of divergent observables 
analyzed in the same theoretical framework provides a 
stringent test of the effective interactions and permits to 
asses the mutual complementarity of reaction and struc- 
ture data for the understanding of exotic configurations 
and decays in those weakly bound nuclei. 

In the SMEC formalism, the subspaces of (quasi-) 
bound (the Q subspace) and scattering (the P subspace) 
states are separated using the projection operator tech- 
nique Q. P subspace contains asymptotic channels, 
made of (N — l)-particle localized states and one nucleon 
in the scattering state. Q subspace contains many-body 
localized states which are build up by both the bound 
state single-particle (s.p.) wave functions and the s.p. 
resonance wave functions. The wave functions in Q and 
P are properly renormalizcd in order to ensure the or- 
thogonality of wave functions in both subspaces. The 
details of the approach can be found in [|||. 

The salient feature of SMEC is that the (quasi-) bound 
many-body states in Q are given by the multiconfigura- 
tional Shell Model (SM) with the realistic effective inter- 
action, providing the internal mixing of configurations. 
The coupling between bound and scattering states is de- 
scribed by the density dependent interaction (DDSM1) 
|5],Q . This interaction provides an external mixing of con- 
figurations via the virtual excitations of particles to the 



continuum states. A subtle balance of external and in- 
ternal configuration mixing explains energies and widths 
of levels, (p,p ) excitation functions, radiative capture 
processes, etc.. 

To generate radial s.p. wave functions in Q and the 
scattering wave functions in P, as a first guess, we use 
the potential of Woods-Saxon (WS) type with the spin- 
orbit : VsoA lr (21 ■ s)r~ 1 df(r)/dr, and Coulomb parts in- 
cluded. A w = 2fm 2 is the pion Compton wavelength 
and f(r) is the spherically symmetrical WS formfac- 
tor. The Coulomb potential Vc is calculated for a uni- 
formly charged sphere. This 'first guess' potential U(r), 
is then self-consistently modified by the residual inter- 
action. The iterative procedure yields the self-consistent 
potential f/( sc ) (r), which depends on channel angular mo- 
mentum and differs significantly from the initial potential 
in the interior of the potential well. For weakly bound 
configurations, also the surface properties of initial po- 
tential are changed by the residual coupling. As a result 
of this self-consistent coupling, the initial SM Hamilto- 
nian in Q becomes energy dependent and its eigenvalues 
become complex for decaying many-body states. 



II. THE FIRST-FORBIDDEN MIRROR p 
DECAYS IN A=17 NUCLEI 

First-forbidden f3 + decay rate (f + ) from the ground 
state (g.s.) J ir = 1/2^ of 17 Ne to the weakly bound state 
J n = 1/2 j in 17 F exhibits an abnormal asymmetry with 
respect to its mirror (3~ decay rate (/ _ ) of 17 N into a 
well bound excited state of 17 ||. This asymmetry has 
been explained either by large asymmetry of radial sizes 
of lsi/2 s.p. orbits involved in bound states of 17 F/ 17 
and 17 Ne/ 17 N J^] or by the charge-dependent effects lead- 
ing to different amplitudes of Tr(0p^ 2 ls1^ 2 )v(0p\^ 2 ) and 
u(Qpy 2 lSy 2 )7T(0py 2 ) components in the g.s. wave func- 
tions of 17 Ne and 17 N, respectively These two quali- 
tative analysis of the mirror asymmetry lack a consistent 
treatment of radial properties of the wave functions which 
are involved in these transitions. 
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FIG. 1. The astrophysical 5*-factor for the reaction 
16 0(p,7) 17 F leading to the states 5/2+ and l/2t is plotted 
as a function of the center of mass energy Eom ¥§■ The ex- 
perimental data are from M. The diffusenes parameter of 
the initial WS potential U(r) is a = 0.55 fm. The depth of 
U (r) is chosen in such a way that the self-consistent potential 
yields the binding energies of proton s.p. orbits 0d 5 /2 and 
lsi/2 at the experimental binding energies of 5/2^ and 1/2+ 
many-body states in 17 F (from MJ). 

The stringent constraint on the diffuseness of the 
mean-field generating radial s.p. wave functions is pro- 
vided by the proton radiative capture cross section (the 
astrophysical S— factor) 16 0(p, 7) 17 F to the 'proton halo 
state' 1/2+, which in the formalism of SMEC can be 
calculated using not only the same wave functions in Q 
and P, but also the same SM interaction, residual cou- 
pling or the initial average potential, as those used in 
the calculation of the decay rate 17 Ne(/3+) 17 F and the 
spectrum of 17 F. Fig. 1 compares measured astrophys- 
ical S— factor for the reaction 16 0(p, 7) 17 F with the re- 
sults of the SMEC calculations using ZBM interaction in 
(Pi/2^5/2 s i/2) shells. The capture cross-section data can 
be well reproduced using U (r) with the surface diffuse- 
ness in the range : a — 0.55 ± 0.05 fm. 

Having determined the surface features of the mean- 
field, we can estimate the charge-dependent effects in 
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dominant contribution to the first- 



forbidden mirror (3 transitions systems. The calcula- 
tion of the nuclear matrix elements and, hence, the first- 
forbidden (3— decay rates are done using the formalism of 
Towner and Hardy Following the analysis of War- 
burton et al. for A ~ 16, the matrix element £ v of 
the time-like piece of the axial current is multiplied by a 
constant factor 1.61 to account for an enhancement due 
to meson-exchange currents. 

The experimental rate /+ for 17 Ne(/3+) 17 F is repro- 
duced by the SMEC calculation with ZBM-F interaction. 
This hybrid interaction reproduces exactly the experi- 
mental energy splitting between the g.s. 5/2+ and the 
first excited state 1/2+ in 17 F. However, for the same in- 
teraction the measured rate /~ for 17 N(/3~) 17 is over- 
predicted by a factor ~ 3. Since the radial dependences 
which are consistent with the proton capture data give an 
excellent fit of both the f3 + decay rate and the spectrum 
of 17 F, the discrepancy for /~ and /+/ /~ is uniquely due 
to the deficiency of ZBM-F interaction to reproduce the 
configuration mixing in 17 and 17 N, i.e. by the charge- 
dependent effects in the Q space SM Hamiltonian. 

How large are these effects? An essential parameter 
here is the amplitude of a component (lsf^ 2 0p^^ 2 ) in the 
g.s. wave functions of 17 Ne and 17 N. Experimental val- 
ues for /~ and / + / '/~ are reproduced if the amplitude 
of (ls^ 2 0p7, 2 ) in g.s. of 17 N is reduced by ~ 30%. In a 
simplest way, this can be achieved by reducing the sep- 
aration of 0d 5 /2 and lsi/2 shells in ZBM-F interaction. 
This new interaction, called ZBM-O*, reproduces also 
well the experimental spacing of 5/2+ and 1/2+ levels in 
17 0. One should notice that this important reduction of 
the (ls^OpT^) amplitude concerns a very small com- 
ponent of the g.s. wave function of 17 Ne and 17 N. The 
dominant component, which is (0d§, 2 0p^ 2 ), changes by 
less than 5% going from ZBM-F interaction to ZBM-O* 
interaction. Also the dominant lp-Oh configuration in 
g.s. 5/2+ and first excited 1/2+ states of 17 F (ZBM-F) 
and 17 (ZBM-O*) is only slightly modified (~ 12%). 
Of course, one should be aware that these estimates 
of charge-dependent effects in SM Hamiltonian can be 
somewhat affected by the chosen ZBM valence space. In 
particular, the absence of OP3/2 and 0d 3 / 2 subshells leads 
to an amplification of the sensitivity in the lsi/2 ~~ > 0pi/2 
contribution to the charge-dependent effects ||. Never- 
theless, for the first time the quantitative extraction of 
charge-dependent effects, separately on the radial prop- 
erties of wave functions and on the configuration mixing 
in mirror systems, became possible. 



III. BINDING ENERGIES IN NEUTRON-RICH 
OXYGEN ISOTOPES 

As a second example of application of the SMEC for- 
malism, we shall consider the correction to the SM g.s. 
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energy due to the continuum coupling. In this section, 



s.p. state. 




FIG. 2. The SMEC energy correction to the g.s. masses 
of oxygen isotopes is calculated in the full sd shell using a 
Widenthal effective interaction llOl. 
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we use the valence space of full sd shell and the ef- 
fective interaction of Wildenthal JhJ. One- and two- 
neutron separation energies (S n and S^n, respectively) 
have been calculated for this isotope chain using SM in 
the same effective space and for the same effective in- 
teraction p"lj] . Fig. 2 shows the above mentioned en- 
ergy correction E corr to the SM masses in the chain of 
neutron-rich isotopes of oxygen. E corr strongly increases 
while approaching the neutron drip line. The 'odd-even 
staggering' reflects a sensitive dependence of this correc- 
tion on the one- neutron separation energy S n , which ex- 
hibits a similar dependence. This alone is however not 
sufficient. The second important effect is related to odd- 
even variation in the density of many-body states in the 
N — 1 system determining the overall number of channels 
by which system N couples to the scattering continuum. 

Particularly large coupling matrix elements correspond 
to the isoscalar couplings between proton and neutron 
fluids. In the example shown in Fig. 2, these couplings 
are absent (number of protons in sd shell equals zero), 
but one expects that the SMEC energy correction to 
SM masses will strongly increase when both fluids are 
present. This may indicate a large shift of the position 
of the neutron drip line between oxygen and fluor. 

SMEC correction to SM masses depends strongly on 
the radial wave functions in Q. An example of this kind 
can be seen in Fig. 2. For oxygen isotopes, the standard 
isospin dependence of the depth of the central potential 
yields Od^/2 s.p. state for neutrons unbound for N < 18. 
With increasing N, the energy of Od^/2 resonance is going 
down and for N — 18 it becomes bound by ~ 200 keV. 
In view of the uncertainty concerning the isospin depen- 
dence of the central potential, we have calculated E corr 
for 26 (N = 18) by keeping 0<i 3 / 2 in the continuum at 
~ 50 keV. The full point for N — 18 gives an idea of the 
sensitivity of E corr to the asymptotic properties of 0d 3 / 2 



FIG. 3. Excitation energy of 2f excited state in even-iV 
isotopes of O calculated in SMEC (solid line) and in SM (the 
dashed line). 
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FIG. 4. Lowest energy states in odd-iV oxygen isotopes 
around 24 O. The dotted line shows the position of one-neutron 
emission threshold. SMEC and SM results are shown with the 
solid and dashed lines, respectively. 

When approaching neutron drip line, the two neutron 
separation energy S^n decreases fast and may become 
smaller than S n . In this case, the dominant continuum 
coupling corresponds to a virtual excitation of nucleon 
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pairs, which are not yet included in SMEC. For that rea- 
son, the SMEC correction close to neutron drip line has 
to be considered as a lower limit of the continuum influ- 
ence on SM masses. 

Independent information about the configuration de- 
pendence of the continuum coupling is provided by the 
position of the 2* state in even-iV oxygen isotopes. In 
Fig. 3, we compare results of SM and SMEC, which are 
obtained using the same model space and interaction in 
Q. One can see that the relative shift of 2^ with respect 
to the g.s. 0^ depends on N, and in some cases the con- 
tinuum coupling shifts this state upwards with respect 
to the SM prediction. However particularly interesting 
is 24 0, for which the dominant g.s. configuration does 
not contain the contribution from 0d 3 / 2 s.p. resonance, 
whereas this resonance is strongly populated in 2f state. 
We can see that the relative shift of 2\ state in SMEC 
calculations for 24 O is particularly strong, reflecting dif- 
ferent radial asymptotic properties of the dominant con- 
figuration. 

A spectacular example of the dependence of the SMEC 
energy correction on the s.p. radial wave functions in Q 
is shown in Fig. 4. Here the lowest energy states of 23 
and 25 O are shown for SMEC (the solid lines) and SM 
(the dashed lines) . Zero of the energy scale corresponds 
to the g.s. in both cases. G.s. of 23 does not contain 
any significant component involving the occupation of 
the s.p. resonans 0d 3 / 2 , and the SMEC correction for 
excited states 5/2+ and 3/2 + close to the threshold (the 
dotted line) are larger than in the g.s. In 25 0, on the 
contrary, 0c? 3 / 2 s.p. resonans is strongly occupied in the 
g.s. and , therefore, excited states l/2 + and 5/2 + in 
SMEC calculations are shifted upwards with respect to 
their position in the SM. 

Not only the absolute energy corrections are large for 
nuclei near the neutron drip-line, as shown in Fig. 2 for 
the g.s. of even-even nuclei, but also the relative shifts 
of excited states with respect to the g.s. can be large 
and of different signs depending on the radial features 
of s.p. orbits involved, the internal configuration mixing 
provided by the SM effective interaction, and the external 
configuration mixing which depends on the size of the 
coupling matrix elements for different channels in [(N — 
l)(gm] J . Systematic experimental investigation of these 
salient tendencies of the multiconfigurational continuum 
shell model remains a challenge for future studies. 



allows to hope for a new wave of exciting future develop- 
ments in the field of exotic nuclei. 
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IV. OUTLOOK 



With the SMEC, the new paradigm is born for the 
microscopic understanding of low-lying properties of un- 
stable nuclei. Fruitful exchanges with the standard SM, 
which nowadays experience a revival in medium-heavy 
nuclei and allows to perform no-core calculations till 16 O, 
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